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We consider the the effect of the electron cloud about a vorton in the CFL-JC" high-density phase 
by numerically solving the ultrarelativistic Thomas-Fermi equation about a toroidal charge. Includ- 
ing electrons removes the electric monopole contribution to the energy, and noticeably decreases the 
equilibrium radius of these stable vortex loops. 



I. INTRODUCTION 

The behavior of matter at very high densities, like 
those found at the core of neutron stars, remains largely 
an open problem. This is due to a lack of reliable meth- 
ods to extract information directly from QCD. But since 
the temperatures relevant for the physics of neutron stars 
is very low (as compared to the Fermi energy) a num- 
ber of observables of interest depend only on properties 
of the low lying excitations of the ground state. In the 
case of dense matter it is believed that those excitations 
are pseudo-Nambu-Goldstone (NG) bosons arising from 
the spontaneous breaking of certain approximate symme- 
tries. The study of these quasiparticles using the meth- 
ods of effective field theory opens up the possibility of a 
more reliable understanding of dense matter. 

The breaking of symmetries in high-density QCD is 
related to quark pairing. At high density, evidence now 
points towards color superconductivity where quarks pair 
in different combinations depending on the temperature 
and chemical potential[l-4]. This pairing uses the BCS 
mechanism, just like electron pairing in conventional elec- 
trical superconductors. However, color superconductiv- 
ity is a much richer phenomenon than electrical super- 
conductivity because quarks carry flavor and color while 
electrons do not — different color superconducting phases 
arise depending on which quarks pair. The nature of the 
NG bosons depends crucially on the pattern of symmetry 
breaking and, consequently, on the kinds of quark pairing 
found in the ground state. 

At asymptotically high density, three-flavor QCD is be- 
lieved to pair in a pattern known as color-flavor-locking 
(CFL), as it pairs all quarks[5]. At lower densities, where 
the QCD coupling is not small, the problem is more com- 
plicated and a plethora of phases have been suggested (for 
a review see [6, 7]). Perhaps the most likely phase at very 
high density is the CFL-if° phase, which seems guaran- 
teed to be the energetically favored phase for some range 
of large densities [8, 9]. This range may include densi- 
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ties relevant for compact astrophysical objects and thus 
may be phenomenologically relevant for extracting pre- 
dictions about extreme environments, like the inner cores 
of neutron stars. 

The NG modes of the CFL phase form an SU(3) octet, 
and are named in direct analogy with the in vacuo SU(3) 
flavor octet, but at high density the "kaons" are the light- 
est particles [10-14]. As they live in an SU(3) octet, their 
behavior is described by a theory resembling familiar chi- 
ral perturbation theory. As the chemical potential ^ is 
lowered and the quark masses become relevant some of 
the NG modes become lighter. When /j, is low enough, the 
neutral "kaons" would have a negative mass-squared, and 
instead condense with some vacuum expectation value, 
turning the CFL phase into CFL-iv:°[8, 9]. This neutral 
kaon condensate provides a superfluid background that is 
the stage for interesting topological defects. The reason 
the neutral, as opposed to the charged, kaons condense is 
the small differences in their mass induced by isospin vi- 
olating quark masses and electromagnetic effects. These 
differences are small so a condensate of charged kaons 
is almost degenerate with the favored neutral kaon con- 
densate. This peculiar situation provides a physical re- 
alization of the so-called "superconducting string" solu- 
tions, suggested in more abstract contexts [15] and also 
appearing in some models describing beyond the Stan- 
dard Model physics[16-19] . Regular global vortices ap- 
pear due to the if" condensation. At their core, the value 
of the condensate vanishes, as it is the case in all vor- 
tices. Since there is a close competition between neutral 
and charged kaon condensates, the charged kaons will 
then condense in the core of the vortex. The con- 
densation of a charged particle leads to electrical super- 
conductivity (not color superconductivity) along the vor- 
tex. Loops of superconducting vortices ( "vortons" ) can 
be stabilized by the supercurrent running around them. 
The mechanism, roughly speaking, is that the energy of 
a vorton has two main contributions. One is the tension 
along their length, which creates an energy that scales 
with its radius R and drives a vorton to be small. How- 
ever, the energy of a vorton includes terms that scale like 
,P /Q^R where J is its angular momentum and Q its elec- 
tric charge. These terms prevent a vorton from collaps- 
ing. This mechanism was invoked in [20] in the context 
of quark dense matter. It was subsequently pointed out 
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[21, 22] that, just like in the case of cosmological vortons, 
the currents necessary to stabihze a vorton could be high 
enough that they would quench the superconductivity of 
the vortex core, destroying its supercurrent and render- 
ing the stability mechanism mute. The quenching effect 
can, however, be counteracted by electric charge. A vor- 
ton with an overall electric charge, in addition, would be 
more stable on account of Coulomb repulsion between 
opposong sides of the loop. 

A more quantitative theory of quark matter vortons 
appeared only recently. In Ref. [23], the equilibrium size 
of vortons containing charge as well current was esti- 
mated by making a variety of approximations and as- 
sumptions. It was found that stabilizing a vorton so that 
its equilibrium radius was several times the thickness of 
its core d and thus could really be identified as a torus, 
required large electric charges {Z ^ 10'*). With such a 
large charge, it is possible that ambient electrons which 
are in the bulk to ensure charge neutrality could be pulled 
close to the vorton, shielding its electric fields, and signif- 
icantly reducing the effectiveness with which the electric 
charge aided in stabilizing the vorton in the first place. 
We call these vortons with neutralizing orbiting electrons 
vortonium to easily distinguish them from unshielded or 
bare vortons. 

It is the purpose of this paper to account for these 
electrons and to find the equilibrium radius of a vorto- 
nium "nucleus" . We use the Thomas-Fermi approxima- 
tion to find the electrons' semiclassical charge distribu- 
tion, allowing them to shield the electric field. Second, 
having found the electron distribution, we can compute 
the energy in the electric field numerically, without mak- 
ing geometrical or multipole approximations. Intuitively, 
when including electrons we should find that vortons have 
a smaller equilibrium radius, because there are exactly 
enough electrons to kill the monopole field. We find that 
making these improvements decreases the equilibrium ra- 
dius of a vorton of given angular momentum and electric 
charge significantly. 
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breaking the approximate U{1)a symmetry that Dirac 
fermions enjoy at asymptotically high density. 

Simple counting suggests that there should be 18 
Nambu-Goldstone mesons, but eight of these are eaten 
by the gluons, which gain masses ^ g^/i, indicating that 
the phase is color superconductive[10]. The condensate 
breaks the electromagnetic U{1)em C SU{3) l x SU (3)^?, 
but there remains an unbroken gauge symmetry that 
arises from a linear combination of the eighth gluon and 
the photon. This unbroken gauge symmetry, which we 
call U{1)q is "mostly" electromagnetism: when Uem is 
small the contribution from the eighth gluon is small, 
the differences between aq and uem are small, and the 
dielectric constant €cfl differs from evacuum by only 
a small amount. All of these corrections are at least 
O {uem) suppressed, and we henceforth refer to the re- 
maining U{\)q symmetry as electromagnetism for con- 
venience. 

After the gluons eat eight of them, ten NG modes re- 
main. Two describe bulk superfluidity that results from 
breaking the baryon and axial symmetries. The other 
eight are an octet of the SU{?>)c+l+r- Because these 
eight "mesons" share group structure with the in vacuo 
mesons that are described by flavor SU(3), these modes 
are named in an analogous manner. The theory contains 
a chiral field E which can take on a vacuum expectation 
value So- The eight mesons tt" are small fluctuations 
about So, so that E = exp (itt^A"//) Eq, where A" are 
the usual 8 Gell-Mann matrices with Tr (A^A^) 2(5*^ 
and / plays a role analogous to the pion decay constant. 

To lowest-order, the effective theory describing the 
CFL phase is described by the Lagrangian [8, 10-13] 
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To understand vortons, we first must understand the 
phase in which they appear. When the quark masses are 
negligible compared to the chemical potential /i, the CFL 
phase is described by the quarks pairing according to [24] 
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where q are the quarks which carry color a, 6 e {1, 2, 3}, 
flavor i,2 G \u^d^s\, and helicity i?, L indices. The di- 
mensionful scale A is known as the gap, and gs is the 
strong coupling constant. This condensate breaks the full 
S'f7(3)c X SU{i)L X SU{-S)r X U{1)b x U{1)a symmetry 
down to SU{3)c+L+R. x Z2 x Z2, where the first remain- 
ing Z2 comes from breaking U{1)b and the second from 



M = diag (m„, m^, is the quark mass matrix, Q — 
|diag (2, —1, —1) is the matrix of the quark charges under 
U{1)q, is the U{\)q gauge field, F its corresponding 
field strength tensor, A, /, and v low-energy constants, 
and pf the Fermi momentum, which we take to be equal 
to the chemical potential 

The leading-order mass term differs from the usual chi- 
ral Lagrangian, which has a term linear in M. A heuristic 
way to understand this difference is that the CFL con- 
densate leaves Z2 remnants of the left- and right-S'C/(3) 
symmetries, so that left- and right-handed quark number 
should be preserved modulo 2, while terms linear in M 
mix left- and right-handedness and violate this conserva- 
tion and are thus forbidden. The mass term included 
in this Lagrangian is not the only O (M^) term, but 
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the other possible terms are higher-order in the power- 
counting scheme developed in Ref. [8]. 

Interestingly, the low-energy constants can be matched 
to perturbative QCD calculations at asymptotic fi. Such 
matching calculations give [10, 12] 
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and A 
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A2. (4) 



Armed with this information, one can easily find the 
masses of the CFL mesons at asymptotically large fi. Ex- 
panding about So = 1, one finds that the O (WP) term 
gives masses [8] 




2 {md + TO„)ms, 



2 mdirris + m„). 



^2 m„(rns + mrf). 
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where iiij = (m| — m|)/(2/i) is the fictitious chemical 
potential which arises from the M commutator terms 
in the time-like derivatives. As /x comes down from 
infinity, the kaons and charged pions get a fictitious 
strangeness chemical potential from the MM* /2pF com- 
mutator terms in the time-like derivatives, and are incen- 
tivized to condense. The other mesons have even larger 
masses, but as they lie on the diagonal of E when the 
other fields vanish, their masses do not get contributions 
from the commutator term and are not encouraged to 
condense. 

Since the neutral kaons are the lightest, they condense 
first. The charged kaons are have almost the same mass 
as the neutral kaons — the difference is dependent on the 
mass difference md — rriu and also on electromagnetic 
effects. However, the interactions between mesons in 
the effective theory strongly disfavors two meson species 
from condensing at once and, if K'^ condensed, electrons 
would have to be sprinkled throughout the bulk with an 
additional cost in energy due to the electrons' Fermi en- 
ergy. 

With all of the other NG modes vanishing, the expec- 
tation value of the neutral kaons is [8] 



cos 



(6) 



where mi^o is given in (5). Since cosine is restricted to be 
in the interval [—1, 1] we see that if rn^i^a/ IJ-ld > ^^ there 
is no neutral kaon condensate at all. From the form of 
the chiral field S = exp (ztt^A"//) Sq, it is clear that the 
value of the kaon expectation value is a compact variable 
whose value only matters modulo 27r, and thus it is not 
surprising that cosine shows up in the expression for the 
K'^ condensate. 



III. BARE VORTONS IN C¥l.-K° 

The background condensate of neutral kaons provides a 
standard superfluid in which we can study global vortex 
strings. As in any standard superfluid, the condensate 
must vanish at a vortex's center in order to remain single- 
valued. One feature of the CFL-/C° superfluid is that 
the other mesons could condense but interactions prevent 
more than one species from condensing. Since the center 
of a superfluid vortex is evacuated, it is natural to expect 
that the cores of these vortices may contain a condensate 
of the next lightest species — the charged kaon. Thus, the 
superfluid vortices have at their core a superconducting 
condensate which may carry electrical current or hold an 
electrical charge. A vorton is a vortex loop with this 
superconducting core. 

A detailed analysis of the stability of a vorton without 
electrons is given in Ref. [23]. There, an in-depth discus- 
sion about the stability of vortons with a flxed electrical 
charge and angular momentum is provided. A vortex nat- 
urally has an associated tension, so a vortex loop would 
naturally shrink to nothing and destroy itself very quickly 
unless stabilized in some manner. When the inner con- 
densate does not break a gauge symmetry, the stability 
of a vorton is easy to understand. Requiring that the 
core contain some flxed number of particles which have a 
preferred density can provide a pressure along the length 
of the vortex. Alternatively, the condensate can have a 
phase which wraps around the length of the vortex. Be- 
cause the condensate must be single valued, this phase 
must be quantized, and is inherently topological. The 
topological winding number N enters quadratically into 
the tension F± as (N/R)'^, while the usual terms in the 
tension are a constant c independent of R (the energy cost 
of suppressing the background condensate minus the sav- 
ings provided by allowing the core to condense). Then 
the energy F ^ 2-kRF± ~ N'^/R + cR has a preferred 
equilibrium radius Rq. 

However, the gauged case is a bit more tricky. It was 
also noted in [23] that due to gauge invariance, the topo- 
logical winding number TV of that wraps around the 
length of the vorton and the magnetic fluxoid threaded 
through a vorton $b, while separate observables, always 
enter the Lagrangian and the energy functional F in just 
the right combination so as to give the angular momen- 
tum J, so that we need not specify them separately but 
instead need only specify J. In field theory, angular mo- 
mentum depends on one space and one time derivative, 
while the charge depends on a time derivative. From the 
form of the derivatives (2) one concludes that the time 
and space derivatives are traded in the energy for Q and 
J/QR. 

In order to get an analytical feel for the energy func- 
tional of a vorton without surrounding electrons, 
Ref. [23] makes several approximations. The first is that 
instead of computing the F directly for a vortex loop, 
instead the tension T± for a straight superconducting 
vortex is used and simply multiplied by the length of 



4 



the vortex 2'kR, which directly neglects curvature effects. 
This contribution to F accounts for all of the space in 
a torus whose radius and thickness are both R. To that 
contribution, the electromagnetic energy is approximated 
by its leading multipole contributions-the electric field 
by the monopole contribution and the magnetic field by 
the dipole contribution. These far-field energies account 
for the space outside a circumscribed sphere of radius 
2R around the inner torus. The apple-core-shaped space 
outside of the torus but inside of the sphere was left un- 
accounted for. Artificially partitioning space and only 
taking only the leading multipole moments both reduce 
the accuracy of F. 

The leading far-field electric part of F is easily com- 
puted and is given by 

The electric part of F that is outside of the vortex core 
but within the torus is given by 

(8) 

In order improve upon the results in [23], we will mini- 
mize a modified F . We will minimize 

F — -Fold — FE-fai — FE-ncai + FtF (9) 

where fold is the main analytical result of [23] and Ftp is 
the energy in the electric fields that we find numerically 
by solving the Thomas-Fermi equations for a toroidal 
charge of radius R, thickness i5, and charge eZ. In the 
next section, we will use the Thomas-Fermi approxima- 
tion to find the semiclassical shape of the electron cloud 
surrounding a toroidal charge. Replacing our approxima- 
tions with Ftp accounts for electrons, avoids partitioning 
space artificially for the electric energy, and avoids the 
multipole approximation for the electric energy. We will 
not change the expressions for the angular momentum or 
other expressions which depend on the electric field. For 
details on the non-Coulomb part of the vorton energy we 
point the reader to Ref. ([23]). 

IV. SHIELDING BY ELECTRONS 

In this section we give the method for numerically cal- 
culating the energy in the electric fields Ft_f, accounting 
for electrons. The Thomas-Fermi approach approximates 
the electron density as 

n = ^p%. (10) 

where pp is a local fermi momentum and we work in units 
where h = c = cq = 1. Gauss' law, 

V^0 = en (11) 



tells us how to compute the electric potential (f) given a 
fixed distribution of charges n. We also include the equi- 
librium condition that the marginal energy cost of mov- 
ing an electron from one point to another should zero, 

d(pF -£(/))= (12) 

where we have used the ultrarelativistic approximation 
\/ p"^ + 'm? ^ pp, as we will find that the electrons' fermi 
momentum is significantly higher than their mass. This 
approximation may not be justified for electrons near the 
edge of the electron cloud, but we will find that most elec- 
trons live near the vorton, and thus are contained in a vol- 
ume O (lO^fm'^), meaning that (for Z ^ 10^) they have 
a Fermi momentum ~ 300 MeV ^ rrif, = .511 MeV and 
the ultrarelativistic approximation is appropriate. We 
then know that 

pp -ecf)^ -ecj)^ (13) 

where 0oo is a constant that describes the electric poten- 
tial at infinity. 

Resolving (10), (11), and (13) for an equation depend- 
ing only on </) yields 

VV=^a2(0-0oo)' (14) 

Since is a constant, V'^(/)oo is zero and may be sub- 
tracted from the left-hand side, and we can absorb <j)ao 
into by a simple change of variables. Equivalently, 
we can pick — 0. This is a manifestation of gauge 
invariance, and picking 0oo = is a choice of gauge. 
Eq. (14) is the ultrarelativistic generalization of the usual 
Thomas-Fermi equation and its solution, with appropri- 
ate boundary conditions, is the central theme of this pa- 
per. The Thomas-Fermi approximation neglects the in- 
teraction among the electrons but, at the high electron 
densities occurring in most of the electron cloud this is 
a very small effect. It also neglects the energy associ- 
ated to gradients of the density. Again, these gradients 
are negligible compared to the momenta of the electrons 
themselves. By solving this nonlinear differential equa- 
tion with = subject to boundary conditions that 
account for the toroidal charge representing the vorton 
and the total charge neutrality of the whole system we 
find a (/> that self-consistently dictates and accounts for 
the distribution of electrons. 

To solve take advantage of the geometry of the situa- 
tion, we solved (14) in toroidal coordinates , defined in 
terms of the usual Cartesian coordinates by 

u — —2 Im arccoth {^{\/ x'^ + ^ iz)/aj 

u = 2 Re arccoth ((v/^^Ty2+iz)/a) (15) 
ip — a,Tctan{y / x) 

where a is the radius of the reference circle (which does 
not exactly correspond to the radius of the vorton) and 
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Re and Im give the real and imaginary parts, respectively. 
We use the fundamental domain 



TT < U < TT 



< V < oo 



TT < ip < TT, (16) 



where a large v corresponds to being near the reference 
circle, and v near corresponds to a very large circle 
that comes close to the z-axis. The surface of the vorton 
occurs at v^ax ~ sinh^^ (R/S)'^ — 1. Our expression 
for the energy of a vorton already accounts for the en- 
ergy in the superconducting core v > Vmax , so we restrict 
to the domain < u < Vmax- Excluding electrons from 
V > Vmax creates an artificial pressure which can make S 
artificially smaller. Numerically we find that this effect 
makes changes to S that are ~10%. Since our problem 
has cylindrical symmetry, the ground-state electron dis- 
tribution will not depend on the axial coordinate ip and 
we suppress it henceforth. 

The Laplacian for these u, v toroidal coordinates is 

1 



This condition matches the electric field to the electro- 
magnetic gauge field inside the vorton, whose form are 
found explicitly in [23]. All the boundaries are subject 
to Neumann conditions (as expected by gauge symme- 
try), so to numerically stabilize the solution we specify 
the gauge (?!)(0, 0) = 0oo = 0. 

The important question now arises: how does the sys- 
tem know how many electrons to use to shield the vor- 
ton at long distances? That is, how does it come about 
that the long-distance field behavior does not include a 
monopole piece? It isn't clear that we have implemented 
this physically desirable constraint. We will now show 
that with the symmetries we have specified, toroidal co- 
ordinates implement this constraint automatically. 

A straightforward manipulation shows that the spher- 
ical coordinate r is given by 



2 2 

r — a 



cosh V + cos u 
cosh V — cos u 



(22) 



V = ^(cos(m) — cosh(t;))x 

[(cos(u) coth(w) — csch(w))9„0 - 
-f (cos(u) - cosh(t-)) {dl + dl) ( 



and so large corresponds to small +v'^. The multi- 



sin(u)9„0 (17) pole expansion of (f>, 



We discretized u and v onto a 32 x 32 grid, and promoted 
(j) to (/)(t) where t is a fictitious time coordinate. Then, we 
numerically solve the coupled ordinary differential equa- 
tion, 
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(18) 



(j){r) = 0o 
corresponds to 

0(u, w) = (f>oo 
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so that if (/>(i/) reaches a temporal fixed point at some 
late time i/, it must also automatically satisfy (14). The 
function multiplying d4>/dt was picked to alleviate some 
numerical instabilities, but the final solution of (14) is 
independent of this choice as long as d4>/dt vanishes. 

The toroidal coordinate u spans a range —n to n. 
The value u = corresponds to the xy— plane with 

+ y'^ > , while u = n corresponds to that same plane 
but with x'^+y'^ < a^. Since our physical setup has parity 
across this plane, we can reduce the region where we need 
to solve (14) to < u < vr, with the boundary conditions 



du(f>{0,v) = du(f>in,v) = 0. 



(19) 



This reduction of the physical space allows us to get a 
more reliable solution while using a grid of the same size. 

Axial symmetry about the z— axis (which corresponds 
to V = 0) entails that the solution must obey 



dy(l){u,Q) = 0. 



(20) 



By symmetry alone we have specified the conditions on 
three out of the four boundaries of the region. 

The final boundary is at the wire. There, we enforce 
that the electric field is that of a wire of linear charge 
density Ze/2nR and circular cross-section nS'^, 



d^4>{u,Vr,iax) = 2 ^^^^ COth(v„ 



(21) 



(24) 

when both u and v are small. We choose <f>aa to be 
zero as a choice of gauge, and d = is enforced by 
parity and axial symmetries. However, it is not obvi- 
ous that we have enforced Qtotai = 0. Part of the puz- 
zle is that all of the surface at infinity is mapped to a 
single point, (m,w) — (0,0). However, on the bound- 
aries u — and v — we have already recognized that 
parity and axial symmetry require du(t>{0,v) = and 
dy(l){u,0) = respectively, meaning that near (0,0) (j) is 
flat in both directions — exactly what we need to guaran- 
tee Qtotai = 0. Charge neutrality is accidentally achieved 
by our other symmetry conditions in these coordinates, 
which makes implementing such a solution easy. 

Unlike the usual application of the Thomas-Fermi 
equation to atoms with spherical symmetry which de- 
pend only on the charge Ze, a dimensionless number, the 
solution for a toroidal charge does not have a simple scal- 
ing behavior, as there are length scale R, 5 and not just Z 
and the electron mass, and the Laplacian does not scale 
simply as the familiar spherical case does. Thus, for ev- 
ery Z, we solve the differential equation for many R and 
5. It is possible to reduce this three-parameter space to 
a two-parameter space but for pedagogical simplicity we 
simply sample the three-dimensional space. We show one 
set of equipotentials from a numerical solution in Fig. 1. 
One can easily verify that far from the charge the equipo- 
tentials appear spherical, slightly closer the charge the 
equipotentials resemble oblate spheroids, near the charge 
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FIG. 1. (color online) Four equipotentials from a typical 
numerical solution to the Thomas-Fermi equation with a 
toroidal charge. 

there are bialy-shaped equipotentials, and most near the 
toroidal charge the equipotentials seem toroidal. 

As a check on the numerical accuracy of the procedure, 
we integrate up the total electric charge outside of the 
toroidal charge Z and hope to find —Z. That is we ask 

-Z= / 27rdA V20(u,z;) (25) 

outside 

where the 27r comes from the axial integration and dA 
is the uw— area differential, including the appropriate Ja- 
cobian. This expectation was satisfied to within 1.8% 
for every choice of parameters; only 2% of the numeri- 
cal solutions had a disagreement between the left- and 
right-hand sides of more than 1%; only a quarter of the 
solutions disagreed by more than 0.5%. Moreover, with 
a finer and finer grid the fulfilling of (25) is better and 
better. 

From the solution 0(u, v) for a given Z, i?, and 5, we 
compute the energy 

FTF{Z,R,b) = J 2^dA \{-V4>f. (26) 

For ease of comparison with the results in [23] 
we computed Ftf for R G [80, 800]fm in steps 
of 20 fm, 5 S [1, 35]fm with steps of I fm, and 
Z e {250, 500, 750} U ([1250, 50000] with steps of 1250) 
and used these data points to construct an interpolating 
function that worked for any set of parameters within 
those ranges and extrapolated to larger 6. Unlike the 
parts of F which care about the chemical potential and 



other parameters which might change depending on the 
environment, Ftf only cares indirectly through R and 
5. Therefore, once Ftf[Z,R,6) is constructed, it can be 
used for any choice of the chemical potential, the gap, 
etc. 

As a way to check the quality of the electron-free ap- 
proximation, we can examine the ratio of the actual en- 
ergy in the Thomas-Fermi calculation to the approxi- 
mated energy, FTF/{FE-fm- + ^^B-noar)- We find that for 
small radii {R < lOOfm), where the electrons can't do 
much shielding, the approximation overestimates by at 
least a factor of 2. At very small radii R ^ d, the electron- 
free approximation may actually be an underestimate, 
but the charge ceases to have an obviously toroidal shape, 
and as other approximations involved in finding the equi- 
librium radius of a vorton breakdown, is not of physical 
interest. With larger radii, where the electrons can shield 
the Coulomb energy well, the approximation tends to be 
off by roughly a factor of 5, but the ratio seems to lose 
any strong R dependence. This comparison is depicted 
in Fig. 2. 
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FIG. 2. The ratio between the numerically computed energy 
Ftf and the approximation Fg.far -I- -Fe-near as a function of 
the vorton radius, for Z G [1250, 40 000] with steps of 1250 
and 5 G [20,35]fm with steps of 1 fm. The exact identity 
of each curve is irrelevant; the generic behavior gives rough 
insight. 

With Ftf{Z, R, S) in hand, we can pick a {Z, J) pair, 
and simply minimize the energy F given in (9). For 
concreteness and ease of comparing to unshielded vor- 
tons, we use the same example parameters as [23], a gap 
A = 66 MeV, an overall chemical potential /i = 450 MeV, 
and quark masses {mu,md,ms) w (2,5,95) MeV. These 
parameters don't appear in Ftf and thus one calcula- 
tion of Ft f empowers us to examine any relevant density 
and gap. We show R/5 as a function of {Z, J) in Fig. 3, 
and we compare the equilibrium radii of shielded and 
unshielded vortons in Fig. 4. With quick examination of 
Fig. 3 it is evident that the effect of electron shielding 
shrinks the equilibrium radius for values of {Z, J) where 
the unshielded calculation is believable (Z > 6000), and 
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FIG. 3. The equilibrium radius _Ro compared to the string 
thickness 5 of a shielded vorton at equilibrium as a function 
of the vorton charge Ze and angular momentum J. 



(Screened R)/(Unscreened R) 




FIG. 4. The ratio R/5 for a vorton with electron shielding 
compared to R/5 for the same vorton with no electron shield- 
ing as a function of Z and J. The rough nature of the small 
Z and J region is due to sampling a rough grid — a finer grid 
would resolve these features. 



that this effect becomes stronger for larger Z. 

One interesting effect that can be seen in Fig. 4 is that 
shielded vortons are actually larger than their unshielded 
counterparts for small Z and J. This behavior fits with 
the low-i? trend seen in Fig. 2 and should be interpreted 
as the electron degeneracy helping prevent the vorton's 
collapse. However, the results in this regime are not 
trustable, as it is precisely the same region where the 
fictitious pressure that arises from restricting the elec- 
trons to u > Vmax shrinks 5 by 20% to 40%. To under- 
stand this area well, one would have to allow electrons 
into the excluded regime in v, and would need to solve 
the Thomas- Fermi equation (14) simultaneously with the 
equations of motion for the two kaon types, which is be- 
yond the scope of this work. 

More careful examination of Fig. 3 shows that not only 
are shielded vortons smaller, but compared to unshielded 
vortons they grow more slowly as one increases Z. This 
is the expected behavior, as to first order putting a posi- 
tive charge in the vorton is cancelled by a shielding elec- 
tron. If this cancellation were exact, a shielded vorton's 
size would be independent of Z . Only when one includes 
Pauli blocking does a shielded vorton's equilibrium radius 
grow. A related final observation is the increasing success 
of the electric shielding by electrons. As Z gets larger, the 
ratio of shielded equilibrium radius and the unshielded 
equilibrium radius drops. This is also expected, as larger 
Z corresponds to larger R, and a larger R enables elec- 
trons to live in the central hole of the vorton more readily, 
shielding the vorton's charge more effectively. 



V. SUMMARY AND OUTLOOK 

In this paper we consider the effect of electron shielding 
on the equilibrium radius of vortons in the CFL- A''^ phase 
at high density. For parameter values where unshielded 
vortons were reasonably large and had a clear toroidal 
shape with an aspect ratio R/5 of a few {Z > 6000), 
the electron shielding of vortonium, as expected, shrinks 
both the radius and aspect ratio for a given Z and also 
diminishes the response of R/S to a change in Z. 

By allowing that portion to be shielded with electrons 
and applying a relativistic Thomas-Fermi approximation, 
we have accounted for the electron cloud surrounding 
such a vorton and have relieved some of uncontrolled 
approximations concerning the toroidal geometry of vor- 
tons in CFL'K^. We have shown that the energy in the 
electric field considered in Ref. [23] is unrealistically op- 
timistic insofar as it drastically changes the equilibrium 
shape of the vorton as compared to the more realistic 
vortonium. The reduced size of vortonium may be un- 
derstood by realizing that the far-fields considered for 
unshielded vortons is predominantly a monopole piece 
while for vortonium it is a quadrupole. 

Two directions seem obvious for pursuing the physics 
of CFL-X^ vortons further. The first is to try to cure all 
of the issues we have touched on in this work: attempt to 



8 



find, in a toroidal geometry, the actual mean-field wave- 
functions of the neutral and charged kaon condensates 
for a vorton, while simultaneously incorporating electron 
shielding and allowing electrons to penetrate the core 
and not artificially excluding them, and could account 
for curvature effects. Work in this direction will almost 
certainly require intensive numerical calculations. 

The other possible direction is to try to estimate the 
vorton production rate in the early times in the neutron 
star life as it goes from a high temperature disordered 
phase to a low temperature CFL-iiT^ phase. This num- 
ber is important to gauge the vorton influence on neutron 
star phenomenology. Along the same lines, one might 
consider the decay rate of these objects through weak 
effects, to try and estimate how long any created popu- 
lation might persist in these extreme environments. This 



line of inquiry would help prioritize any other work on 
vortons at high density, and would provide the most di- 
rect route to understanding if vortons might contribute 
to observable astrophysical phenomena. 
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